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We present Monte Carlo simulations of a dense polymer melt which shows glass-transition-like 
slowing-down upon cooling, as well as a build up of nematic order. At small wave vectors q this model 
' system shows excess scattering similar to that recently reported for light-scattering experiments on 

some polymeric and molecular glass-forming liquids. For our model system we can provide clear 
^"*^ ■ evidence that this excess scattering is due to the onset of short-range nematic order and not directly 

related to the glass transition. 

<n : 

The glassy freezing of supercooled fluids has been a longstanding topic of research in a wide range of disciplines from 
engineering to theoretical physics. Despite these efforts, no agreement concerning the basic physical mechanisms of the 
transition has emerged . A particularly puzzling experimental finding is the onset of some large-scale correlations 
in various supercooled molecular and polymeric glass- forming fluids . The evidence for theses correlations consists 
+Jj \ of excess scattering at small wave numbers q seen in light scattering experiments, which by far exceeds the value to 
be expected from the compressibility of the system. The g-range corresponds to distances of several hundreds to 
\ thousands of A, depending on the material. 

This finding is completely unexpected from the point of view of theories that treat the glass transition as a purely 
kinetic phenomenon, such as e.g., mode coupling theory [Q. The only relevant length scale in this theory is the typical 
inter-particle distance, that is the diameter of the "cages" , the decay of which is described self-consistently by the 
\ theory. 

There are, however, theories [^| postulating an underlying second-order phase transition (possibly kinetically 
masked), and such theories would imply the existence of a growing - and ultimately diverging - intrinsic length 
scale associated with the glass transition, although most existing theories along such lines seem questionable ||. 
Furthermore, direct evidence for such a length scale from computer simulations of simple models undergoing a glass 
transition is rather scarce jS 1 1 and any such lengths detected so far are quite small and can not be connected with 
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the large lengths seen in scattering experiments [g-Q far above the glass transition. 

In the present work we show that this surprising scattering behavior of some supercooled liquids can possibly be 
| attributed to local orientational ordering in these systems. We will present a Monte Carlo simulation of the bond 
fluctuation model of a polymer melt Jll|Jl^] with some orientational ordering tendency. In this lattice model each 
repeat unit of the chains occupies the eight vertices of a unit cube on the three dimensional simple cubic lattice. 
Each unit is connected to its neighbors along the chain with bonds fluctuating in length between 2 and y/lO. The 
stochastic dynamics - single monomer hopping or slithering snake moves - obeys connectivity and excluded volume 
constraints |l2|,|l3|. Using a Hamiltonian that favors long bonds, this model has been extensively used to study the 
glass transition in polymer melts and has been shown to reproduce much of the phenomenology that is observed in 
I ■ experiments E3. 

For the present study we now employ a Hamiltonian which incorporates both bond length and bond angle energies 
]l5fl . The bond length energy favors short bonds, 6 = 2, and the bond angle energy favors stretched angles, 8 = 180°. 
Details are given in (ll| . This model system is cooled down in a stepwise fashion and is equilibrated at each temperature 
with the highly efficient slithering snake algorithm |l7j . After equilibration the dynamics of the chains is studied with 
the single monomer hopping algorithm which is known to reproduce the Rouse dynamics for short chains in the melt 



To study the freezing behavior we analyzed the mean square displacements of the centers of mass of the chains in 



the long time limit. A measure of the longest relaxation time of the chains can be defined as the time for which this 
mean square displacement is equal to the mean square radius of gyration of the chains in the melt. Figure 1 shows a 
plot of this relaxation time, divided by its infinite temperature limit, for chains of length N = 20 in the form of an 
inverse activation plot. In the temperature range investigated we observe a slowing down of the chain dynamics by 
two orders of magnitude which is described well by a Vogel-Fulcher law with a Vogel-Fulcher temperature of about 
To = 0.144. This result would suggest the existence of a glass transition in the melt somewhere around T = 0.2. 

In Fig. 2 we show the scattering intensity S(q) at small wave vectors. The data are for a system of 13732 chains 
of length N = 10 in a cubic box of linear size L = 130. All curves are obtained by averaging S(q) over about 100 



1 



statistically independent configurations. The curve at T = 0.376 is flat and its value is the one to be expected from the 
compressibility of the model system Excess scattering can be found starting at temperatures around T = 0.313 
and its intensity increases by almost a factor of 4 upon decreasing T to 0.219. The g-range where this effect occurs 
corresponds to distances of 10 — 60 lattice constants. Since one effective bond of our model corresponds to a group of 
about five successive chemical bonds |llj-|l4|], one estimates that the lattice spacing corresponds to about 2 A. Thus 
qualitatively this effect is similar to what was observed experimentally in certain glass- forming materials Q. 

In this temperature range, however, another property of our model system comes into play: the Hamiltonian we used 
favors stretched bond angles and short bonds. The conformations of the polymer chains therefore gradually change 
from the Gaussian coil limit at high temperatures to the rigid rod limit at very low temperatures. In the temperature 
region where they can be considered to be almost rod-like objects they should therefore display a tendency for liquid 
crystalline ordering. According to Onsager theory Jl8| this tendency should depend on the aspect ratio (length/ width) 
of the chains in their rod-like state. Figure 3 shows that we observe exactly this behavior. At temperatures around 
T « 0.27 the system starts to develop nematic order, as shown by the behavior of the nematic order parameter P2, 
which is defined as the largest eigenvalue of the Saupe tensor 

Q a p = -j f~ ( u ia u i0 } - ^8 a0 ] , (1) 

where lit is a unit vector along bond i and the sum is over all bonds in the system. The chains of length N = 10 
show only a small increase in the value of the order parameter which is due to a build-up of merely local nematic 
order |l6| ]. The nonzero value of the order parameter for N = 10 thus is just a finite size effect. As demonstrated 
below, this local ordering is, however, sufficient to produce the excess scattering shown in Fig. 2 for the chains of 
length N — 10. On the other hand the chains of length N — 20 with their larger aspect ratio develop long-range 
nematic liquid crystalline order at temperatures below T — 0.26 (l6| . For this system the excess scattering occurs in 
the same g-range and temperature range as for the one with N — 10 but has a larger amplitude. In this case the 
isotropic-nematic transition preempts the glass transition, which would occur at a lower temperature. 

To demonstrate the connection between the excess scattering and the local nematic order, we first of all note 
that the temperature dependence of the excess scattering is the same as that of the nematic short range order. To 
make this link more quantitative we consider the following approximation to the scattering function. At T = the 
ground state of the system consists of perfectly aligned rigid rods. In this case the scattering function exactly factors 
into a contribution from the positional correlation of the rods, S cm (q), and a contribution from the intra- molecular 
correlations, which is just the form factor of the rods, f{q): 




S{q) = S cm {q)f{q), (2) 



where the form factor of the rods is given by 



Here N scatt is the number of scatterers per rod and I is its length. We now approximate our scattering intensity at 
finite temperatures by Eq. (2), neglecting fluctuations of the local director and the finite flexibility of our chains. To 
calculate the form factor of the rods equivalent to our chains we set I equal to the end-to-end distance of the chains 
and set N scaU equal to \J {R 2 } / (b 2 ), where v/ {b 2 ) is the mean bond length. If we divide the scattering function 
for the systems with chains of length N = 10, presented in Fig. 2, by the form factor of the rods determined in this 
way (Fig. 4), we see that the excess scattering completely vanishes. This clearly shows that the excess scattering is 
due to the occurrence of orientational order in our system and that it is even sufficient to have only local orientational 
order to reproduce this effect. 

Thus we have demonstrated that for our model system which exhibits on the one hand slowing down typical for 
a polymer melt undergoing a thermal glass transition and on the other hand a transition from random coil to semi- 
flexible and finally rigid rod behavior of the chains and an accompanying (local) nematic ordering, the excess scattering 
intensity at small wave vectors can be unambiguously ascribed to the development of the local orientational order 
and is not directly connected with the glass transition in the system. Since it is sufficient to have a tendency for local 
nematic ordering for this effect to occur, the question arises whether this mechanism may also be responsible for the 
excess scattering experimentally observed for some glass-forming systems. All of these systems are either polymers 
or consist of highly anisotropic molecular liquids with the accompanying highly anisotropic interaction. While in our 
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model the local orientational order is associated with local chain stretching and is of uniaxial nematic type, it is very 
plausible that other systems may have other types of local orientational ordering, developing from a dense packing 
of anisotropic molecular groups. This should, however, show up in qualitatively similar excess scattering. It would 
therefore be very interesting to reexamine the local orientational properties of these experimental systems and look 
at the effects of external orienting fields on the excess scattering. 
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FIG. 1. Longest relaxation time of chains of length N — 20 as a function of temperature. The dashed line shows a 
Vogel-Fulcher fit giving a Vogel-Fulcher temperature of T = 0.144. 



4 




FIG. 2. Scattering intensity at small wave vectors for a system of chains of length N = 10. The different curves 
temperature range over which one observes the build-up of local orientational order. 
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FIG. 3. Nematic order parameter P2 as a function of temperature. The open circles denote systems of 6866 chains of length 
N — 20 and the filled circles denote systems of 13732 chains of length N = 10. The linear dimension of the box is L — 130 in 
both cases. 
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FIG. 4. Scattering intensity at small wave vectors divided by an equivalent rigid rod form factor of the chains for different 
temperatures as shown in the figure. The results are for chains of length N — 10. 
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